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Abstract 

Based on kinetic equations for the density matrix, drift-diffusion equations are derived for a 
two-dimensional electron gas with Rashba spin-orbit coupling. Universal results are obtained for 
the weak coupling case. Most interesting is the observation that with increasing spin-orbit coupling 
strengths there is a sharp transition between spin diffusion and ballistic spin transport. For strong 
spin-orbit coupling, when the elastic scattering time is much larger than the spin relaxation time, 
undamped spin-coherent waves are identified. The existence of these long-lived spin-coherent states 
is confirmed by exact analytical results obtained from microscopic kinetic equations valid in the 
ballistic regime. 

PACS numbers: 72.25.-b, 72.10.-d, 72.15.Gd 



1 



I. INTRODUCTION 



One of the major challenges of spintronics is the generation and manipulation of a spin 
polarization by exclusively electronic means in nonmagnetic semiconductors at room temper- 
ature. A lively discussed phenomenon in this field is the spin-Hall effect , the theoretical 
description of which is still very controversial. In the experimental analysis of this effect 
2, 3), |j, 0], an electric-field induced spin accumulation is observed near the edges of the sam- 
ple. A widespread idea to understand this phenomenon relies on the notion of a spin current 



oriented transverse to the applied electric field. 



2J This current is assumed to induce a spin 



accumulation in the spin-Hall experiment. However, this seemingly clear 
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picture 
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became the subject of serious controversies in the literature. 
The main problem results from the fact that the spin is not conserved in spin-orbit coupled 
systems so that the definition of a related spin current becomes ambiguous. Alternative 



approaches seem to be more suitable, which avoid the identification of an appropriate spin 
current. Most profound approaches of this kind start from a microscopic theory for the 



spin-density matrix 
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22j or Keldysh Green functions 



23 



24 



followed by an analysis of its long- wavelength and low- frequency limit. These treatments 
provide a complete physical description of spin-related phenomena including accumulation, 
diffusion, and relaxation of spin as well as magnetoelectric and spin-galvanic effects. [3] 
Based on this approach, it was recently concluded that as soon as spins are polarized, the 
spin-orbit interaction leads to spin precession and related spatial oscillations of the spin ac- 
cumulation. [25|] A specific difficulty of the phenomenological drift-diffusion approach is the 
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271 ] to solve the differential equations 



formulation of appropriate boundary conditions 
for the position and time dependent spin and charge densities. 

Another challenge in the field of spintronics refers to mechanisms that allow long spin 
relaxation times. In many situations of practical relevance, the main spin relaxation mech- 
anism is due to spin-orbit interaction (Dyakonov-Perel spin relaxation), which leads to an 
effective internal magnetic field that causes spin precession in the plane perpendicular to 
the magnetic field. Momentum scattering on impurities randomly reorient the respective 
precession axis so that an averaged spin dephasing results. From an application point of 
view, the suppression of this spin dephasing in spintronic devices is a demand of high pri- 
ority. An enhancement of the spin relaxation time has been predicted to occur in spatial 
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regions near the edges of a spin-polarized stripe. [28] A suppression of spin relaxation is 
also observed in samples with an initial periodic spin profile. [28] The related spin-coherent 
standing wave has a several times longer spin relaxation time than a homogeneous electron 
spin polarization. The phase as well as the long spin relaxation time of these novel spin 
waves could be important to reading, writing, and manipulating information in spintronic 
devices. 

In this paper, we focus on coupled spin-charge drift-diffusion equations, which allow a 
treatment of spin-related phenomena that are interesting both from a theoretical and exper- 
imental point of view. Starting from microscopic transport equations for the spin-density 
matrix of a two-dimensional electron gas (2DEG) with spin-orbit coupling of the Rashba 
type, we derive in an exact manner drift-diffusion equations for coupled spin-charge exci- 
tations. These equations, which are treated both for weak and strong spin-orbit coupling, 
allow a due consideration of both the field-induced spin accumulation and spin-coherent ex- 
citations. In the weak-coupling limit, the basic results are universal, i.e., they qualitatively 
agree with findings recently obtained for completely other systems, e.g., the spin transport 
of small polarons. 29] Whereas the weak-coupling case has received already considerable 



interest in the literature 
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30], there are only few considerations of effects due to 



strong spin-orbit coupling. As shown below completely new and interesting phenomena are 
predicted to appear. We mention the transition from the diffusive to the ballistic spin trans- 
port regime and the occurrence of spin-coherent standing waves. In this paper, emphasis is 
put on the interesting regime of strong spin-orbit interaction. 



II. BASIC THEORY 



Coupled spin and charge excitations are treated by an effective-mass Hamiltonian that 
describes electrons and their spin of a 2DEG with short-range spin-independent elastic scat- 
tering on impurities. In addition, the carriers are coupled to each other via spin-orbit 
interaction of Rashba type. To keep the approach transparent, we postpone the considera- 
tion of effects due to external electric and magnetic fields to a forthcoming publication. The 
second quantized form of the Hamiltonian reads 

H = } u a[ x [e h - e F ] a kX - {fiw k ■ cr AA ,) a\ x a k \' + u ^ a k x a k'X, (1) 

k,X fc,A,A' k,k> A 
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where a) kX (a k x) denote the creation (annihilation) operators with in-plane quasi-momentum 
k = (k x , k y ,0) and spin A. In Eq. Bp denotes the Fermi energy, cr the vector of Pauli 
matrices, and u the strength of the 'white-noise' elastic impurity scattering, which gives rise 
to the momentum relaxation time r calculated from 

1 2>ku 2 



J2 5 ( £ «-£k)- (2) 



r h 

fe' 

Other quantities in Eq. flTJ are defined by 

h 2 k 2 f~i , ma fnS 

e k = -z — , ^fc = — {Kxk), K = —-e z , (3) 
2m m nr 

with a denoting the spin-orbit coupling constant and m the effective mass. 

The complete information about all physical quantities of interest are included in the 
charge and spin components of the spin-density matrix 

ti(k,k'\t) = (a[ x a k , x ,) t , (4) 

which is more conveniently expressed in the following manner 

f(k, K | t) = Y, fx (k, K\t), /(fe, K\t)=Y fx'(k, K | t)<T\\i . (5) 

A A, A' 

The k dependence, introduced by the replacements k — > k + k/2, k' —> k — k/2, refers to a 
possible spatial charge and/or spin inhomogeneity and leads to a coupling between charge- 
and spin-degrees of freedom. Treating elastic impurity scattering in the Born approximation 
and keeping only the lowest-order contributions of the spin-orbit interaction in the collision 
integral, we obtain the following Laplace-transformed kinetic equations [14] 

sf - -(« • k)f - -K ■ (f x «) = -(7 ~f) + fo, (6) 

mm t 

sf + 2(uj k x /) - -(k • k)f + —(K x K )f 

m m 

t r oe k t oe k 

where fo and fo denote the initial charge and spin distribution, respectively. The cross-line 
over fc-dependent quantities denotes an integration over the angle tp of the in-plane vector 
k and s is the variable of the Laplace transformation that replaces the time t. The kinetic 
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Eqs. <Q and (J7J have been obtained under the condition h 2 K 2 /(mEF) <C 1, when the spin 
contribution on the left-hand side of Eq. ([6]) can be considered as a small correction. 

To proceed we derive coupled spin-charge drift-diffusion equations, which are completely 
in accordance with the kinetic Eqs. and ([7j) in the long- wavelength limit. The Eqs. §5§ 
and ([70 are written in an more convenient and explicit form as a set of four coupled linear 
equations 

af + iVL(q x f y - q y f x ) = R 

h df 

af x + 2Qcos((p)f z - iQq y f = R x + Qsin((p)-— 

af y + 2Q sm((f)f z + iVtq x f = R y - Q cos(<p)- — 
af, - 2tl [cos{ip)f x + sm{<p)f y ] = R z . (8) 

Here, we introduced the abbreviations Q = ujj-t, q x y = K XjV /k and 

k 

a = cr - iVL— [q x cos((p) + q y sin(^)] , a = sr + 1. (9) 
K 

The quantities on the right-hand side of Eq. (JSJ) are given by 

R = J+rf , R z = J z + rf 0z , (10) 

fi () faj o 

Rx = f x + rf Qx —Qfsm(ip), R y = f+ T f 0y + -—Qfcos(tp). (11) 

After solving the set of linear Eqs. (|HJ) for the components of the density matrix, the solution 
is expanded up to second order in q x and q y . The final analytical integration over the angle 
ip results in closed equations for the physical quantities of interest namely / and /. The 
charge density /, is expressed by f(£k, q \ s) = n(eu)F(q \ s), where n(efc) denotes the Fermi 
distribution function for carriers in thermal equilibrium. This ansatz is fundamental and 
completely in line with our basic assumption concerning the relevant time scale of drift- 
diffusion effects, which should last long enough to guarantee that carriers quickly reestablish 
the thermal distribution. For these relaxation processes, the energy relaxation time t £ is 
responsible. Usually, the inelastic scattering time is much larger than r so that the time 
scale t > t £ that can be accounted for by the widespread drift-diffusion approach is restricted 
to t ^> r (or sr 1). The calculation outlined above is straightforward and leads to the 
result 

d -nJ^j = ^—(ejy 1 ^nj^td^^—iejy, (12) 



where the prime denotes the derivative with respect to This result is used to derive 
equations for all components of the density matrix, which decouple into two sets of equations 
for the components / and i(n x f) z = f r on the one hand and f z and in ■ f = f d on the 
other hand. We obtain 

[ S +D Q ( S ) K 2 ]J+V(s)J r = f , 



s + 



s + 



T s ±(s) 
1 

1 



T S _L{S) 

with the following diffusion coefficients 

Do(s) = 



+ D r (s)K 2 
■ + D z (s)k- 
+ D d {s)K 2 



f r - V r (s)K 2 f = f r0 , 
Iz - Vz(s)J d = f z0 , 



f d - Vd (s)K 2 f z = f ( 



do- 



D 



(J, 



2 • 



D = v 2 t/2, 



D = (T 6 + 6a 4 ^ 2 + 36affi 4 + A8Q 6 D 



DJs) 



12ft 2 



( T 2 ( ( T 2 +4ft 2 ) 2 ( ( T 2 + 2ft 2 ) 

Uq - 6(x 4 ft 2 + 12a 2 ft 4 + 16ft £ 



DJs) 



( x 2 ( ( T 2 + 4ft 2 ) 2 (a 2 + 2ft 2 ) 



spin relaxation times 



2ft 2 



0-Q 



t sz (s) = r s a /2, l/r s = ADK 2 = 2ft 2 /r 



and coupling terms 



V(s) 



v Ka 2 + 2ft 2 
^¥cr 2 + 4ft 2 ' 



V z {s) = v 



2ft 



a 2 + 4ft 2 ' 



(T k 



VJs) 



1 + 



a 2 + 4ft 2 e k n r 
a 2 + 2ft 2 n 



2tta 2 



(13) 
(14) 

(15) 

(16) 

(17) 

(18) 
(19) 

(20) 

(21) 
(22) 



(a 2 + 4ft 2 ) (a 2 + 2ft 2 )' 

The velocity v is given by hk/m. All spin relaxation times are proportional to 1/r (1/t s = 
2c<j 2 r), which is in accordance with the Dyakonov-Perel spin-relaxation mechanism. The 
basic solution in Eqs. ffl3l) to ffTB"]) allows in principle the consideration of non-Markovian spin- 
related phenomena because all transport coefficients depend on the variable s of the Laplace 
transformation. Moreover, in the strong-coupling limit (ft ^> 1), when the inequality r ^> t s 
is satisfied, the spin diffusion exhibits principal a non-Markovian character. Unfortunately, 
it turns out that the full consideration of the complicated time dependence in the drift- 
diffusion Eqs. f fT3l) to fTTBT) is not an passable route. A discussion of this mathematical 



-fz( k , K ) 



(23) 



subtlety is postponed to a forthcoming paper. In the weak-coupling limit Q <C 1, under the 
condition r<T„ the spin evolution proceeds completely in a Markovian manner. 

A striking feature of our basic result is the decomposition of the transport Eqs. (TT3"j) to 
ffT6"j) into two sets of coupled equations for (i) the charge density / and the transverse spin 
component f r and (ii) for f z and the longitudinal spin component f j. This feature has 
already been identified for spin-orbit coupled extended electronic states [20( and for hopping 
transitions of small polarons. [29] As a consequence of this decomposition, the internal 
effective magnetic field H, which is associated with the spin-orbit interaction, is only due 
to the spin components f z and f d but not to f r . From the Maxwell equation, we obtain 
H(r) = -Vip(r) with 

i a/texp \ — lh x x — tHy'y — is, \ z \\ 

k 

where \xb denotes the Bohr magnaton. The measurable magnetic properties are therefore 
completely determined by Eqs. (j!5p and (1161) . With other words, the carrier diffusion leaves 
the spin-induced magnetization unchanged. The situation changes drastically, when an 
electric field is applied. In this case, all components of the density matrix couple to each 
other 29j so that an inhomogeneity in the distribution of charge carriers may also induce 
a magnetization. Recently this magneto-electric effect has been thoroughly investigated in 
the literature for semiconductors with spin-orbit interaction. 

III. WEAK SPIN-ORBIT COUPLING 

The character of spin transport and spin diffusion strongly depends on the strength of 
the spin-orbit coupling, which is expressed by the dimensionless parameter Q = Kl. In 
this Section, we focus on the weak-coupling case Q <C 1, when the mean free path I is 
much shorter than the characteristic length K~ l of the spin-orbit coupling. Accordingly, 
the drift-diffusion Eqs. (TT31) to (TTBT) are expressed by 

— fiK 

( s + D K 2 )f + —f r = f , (24) 
m 

s+ I + CK A 7r _M(£^ K7=/rt , (25) 
t s J m n 

+ Dk 2 ^J J z - ADKj d = f z0 , (26) 



(s + ^ + Dk 2 ^ f d - ADKK 2 f z = f d0 , 



(27) 



with the spin relaxation time 

T 1 

Ts = 2W = ADR!' (28) 
These equations formally agree with results obtained for the hopping transport of small 



polarons. [29|] The similarity becomes complete, when the Hall mobility uh in the lat- 
ter approach is replaced by the quantity (er s /m)(£fcn)'/n. Therefore, we conclude that 
independent of the transport mechanism, the kinetic Eqs. (12^1) to (12?]) are universal for 
weakly spin-orbit coupled carriers. This conclusion is confirmed by similar studies in the 
literature. Q 

The weak-coupling limit is most important for the hopping regime, where the hopping 
length must be identified with the mean free path. For the transport of small polarons, this 
length is the lattice constant a so that the weak-coupling condition Ka <C 1 is almost always 
satisfied. 

Another conclusion follows directly from Eq. (|28|) . For weak spin-orbit coupling, the 
spin relaxation time r s is much larger than r. Consequently, the charge degree of freedom 
approaches more quickly the state of equilibrium than the spin degree of freedom so that 
the frequency response of spin-related quantities to external perturbations does not exhibit 
a resonance character. According to Eq. (T25]) . an evolution of the transverse spin component 
f r is expected that follows adiabatically the time variation of / when t s < t £ . 

The solution of the set of Eqs. fl24l) to (12"7|) has already been thoroughly studied previ- 



ously. [19j From the coupled Eqs. (I21I) and (I25I) . an interesting spin-related effect has been 



predicted in the literature 20] namely the spontaneous splitting of a density packet into two 



counter propagating packets with opposite spins. Other results, which are obtained from 



the universal Eqs. (I24I) to (|27|) were derived and discussed in Ref. [29]. 

Finally, let us add a remark concerning the component of the spin-density vector /, 
which is proportional to n. To lowest and first-order in k, we obtain for the charge density 
/ = n(sk)/s. Using this result and accounting for Eqs. (12), (fTSjl . and ( fT6|) . we arrive at a 
solution 

(^) - _ 2 ,_ , oo.L- , 7\ ( 29 ) 



J ft 



(4bt + 2Q 2 (2st+ 1) 



a s 

that was derived already in a previous paper. [14( This particular contribution gives rise to 
a dissipationless spin current in the ground state [l4| that reflects the time evolution of the 
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spin accumulation. Due to the pole structure of this expression, resonances are predicted to 
occur at characteristic frequencies. [14|, [30] 



IV. STRONG SPIN-ORBIT COUPLING 

For strong spin-orbit coupling Q 3> 1, the character of the set of Eqs. (fl3|) and (fl4j) for 
/ and f r does not change remarkably. We obtain D r = 3D/2 and t s ±_ = r, and therefore 
Ts_i_ <C t £ so that in the regime t ^> r, which is accessible by the drift-diffusion approach, the 
transverse spin component f r adiabatically follows the evolution of the charge density /. 
To confirm this assertion, we analyze Eq. ([Ml) in the strong-coupling limit. Neglecting the 
small contributions s and 3Dk 2 /2 on the left-hand side of Eq. (TM"|) . we immediately obtain 

7 r (e k , r,t)) = ^(l + *£) p(e k , r, t), (30) 

where f(ek,r,t) denotes the exact solution of Eq. fflBl with V(s) = 0. 

While the set of Eqs. ( TLB"]) and (fl^|) for the quantities / and f r provide expected solutions 
also in the strong coupling limit, we observe a dramatic and abrupt change in the behavior of 
the spin components f z and f d described by Eqs. ( IT5l) and (|T6l) . As seen from Eq. ( fl9l) . the 
diffusion coefficient D z vanishes at = 1/ vT2 so that for this coupling strength the physical 
picture of diffusion completely collapses. With further increasing Q, the diffusion coefficient 
D z becomes negative, which indicates an instability in the spin subsystem implying that the 
applicability of the kinetic Eqs. (fT5l) and ffl6l) becomes problematic. According to Eq. (T28l) . 
we have t s < r in the strong-coupling limit (Q > 1). Consequently, the spin diffusion 
evolves on a time scale that is shorter than the elastic scattering time r so that the kinetic 
Eqs. (1151) and (f!6l) describe in principle a non-Markovian behavior. Consequently, only in 
the steady state (s — > 0), one can strictly speak about a Markovian spin diffusion in this 
regime. As shown below, another peculiarity arises namely that the quantities f z and f d 
rapidly change on a length scale of the order of the mean free path /, a result which does not 
correspond to the n expansion of Eqs. §6§ and ([7]). Putting all together, we conclude that 
with increasing coupling strength Q, a sharp transition appears from a diffusive behavior to 
a ballistic regime. This transition could be due both to an increasing spin-orbit coupling K 
and an increasing relaxation time r. 

Since the principal applicability of the drift-diffusion approach for strongly spin-orbit 
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coupled carriers becomes doubtful, we return to the basic Eqs. <Q and ([7]) and consider its 
exact solutions, which are valid in the pure ballistic regime (r — > oo) that corresponds to 
the extreme strong coupling limit. First, we treat the ballistic evolution of a spin packet 
that is initially (t = 0) created at r = 0. In the limit r — > oo, the kinetic Eqs. and ([7]) 
are solved by 

fzO 



/3(n, s) 



(31) 



P{k, s)+4u 2 k /(3{K, S y 

where v = hk/m denotes the electron velocity and where an additional contribution has been 
omitted that disappears after the integration over the angle <p. The back transformation to 
the r, t representation is easily accomplished and we get 



f z (r,t) = f z0 cos(2Kr)5(r 2 -vH 2 ). 



(32) 



This solution expresses the pure ballistic character of the spin transport. 

The front of the conic spin packet moves with 
the constant velocity v whereas its amplitude os- 
cillates with a frequency 2Kv due to spin rotation. 
This oscillation manifests itself by a spatial peri- 
odic pattern, which in the time domain describes 
spin precession with the Lamor frequency 2o>fc.Let 
us treat another example namely the half plane 
y > with a given spin polarization at the bound- 
ary y = 0. This steady-state boundary condition 
is described by a term f z o/s on the right-hand side 
). Again, we start from the solution in 
Eq. (13TT) , restrict to the y dependence and perform 




20 30 40 



60 



FIG. 1: Spin polarization vf z /f z o as a tunc 
tion of £ = 2Ky calculated for the half space r p 
y > from Eq. (33). q ' 



the inverse Fourier and Laplace transformations. The solution 



Ll 

fzO 



2u k t - £ 
2ttv 



6(2^-0 / d V 



cos [{{2u k t - OV + 1)] 



^((2uo k t - Or] 



(33) 



manifests a spin front moving with velocity v in the y direction (£ = 2Ky). In contrast to the 
point source in Eq. (|32|) . one observes non- vanishing oscillatory spin excitations also in the 
interior. Spin-coherent waves are formed as shown in Fig. [TJ The physical origin of this spin 
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pattern is due to the instability of the system expressed by a negative diffusion coefficient 
(D z < 0). In this regime, spin diffusion has the tendency to strengthen an initial spin 
inhomogeneity, whereas the competition with spin relaxation leads to the spatial oscillatory 
spin pattern. 

We point out that the exact solutions of the kinetic Eqs. (JSl) and (j7j) in the ballistic 
regime describe an oscillation spin pattern, which changes on a length scale K~ x that is 
much smaller than the mean free path I (Q = Kl ^> 1). These exact results derived for the 
ballistic regime will help to recognize the solution of the drift-diffusion Eqs. (1151) and ({TB"]) 
in the limit of strong spin-orbit coupling. 

We start our analysis of these equations by calculating its possible eigenmodes. Treating 
the determinant of Eqs. (|15|) and (fTBj) in the limit s — > 0, we see that there is no oscil- 
latory mode in the weak-coupling regime f2 < l/(3\/3). However, for coupling strengths 
in the interval l/(3y / 3) < < l/(2y / 3), there are two solutions with vanishing imaginary 
parts that describe completely undamped oscillations. For even higher spin-orbit interaction 
l/(2\/3) < O, only one oscillatory solution is obtained. To illustrate this behavior let us 
treat a stripe of width 2L parallel to the x axis (-L < y < L). A permanent spin polar- 
ization should be provided at the boundaries [f z (y = ±L, s) = f z o/s]. For this special case, 
the inverse Fourier transformation of Eqs. (fT5l) and ({TB"]) leads to the following dimensionless 
differential equations 



with F d = f d / K. The derivative refers to £ = 2Ky and the parameters are given by 



d z f s ~ ®zf z + ldF d = 0, 
d d F' d ' - a d F d - j z f z = 0, 



(34) 



(35) 



a z = a sr + 4tt 2 , a d = a^sr + 2fi 2 (2sr+ 1), (36) 




For the boundary conditions f z o ^ and f y o = 0, we obtain 

j f z o K±(a z — d z K%) sinh /t 2 £o cosh k±£ — K 2 (a 2 — d z n\) sinh Ki£ cosh 

s Ki(a z — d z H^) sinh ^2^0 cosh Ki£o — ^2(^2 — d z K%) sinh /ti£o cosh ' 



11 



with £ = 2KL. ki j2 are solutions of the biquadratic equation 



(d z n 2 - a z )(d d K 2 - a d ) + ^ z ^ d ^ = 0, 



(40) 



with the plus sign. The remaining inverse Laplace transformation requires a numerical 
integration. 

First, let us treat the steady-state solution (s = 0) of Eq. (|39|) . 

Figure 2 shows the spin polarization f z 
function of £. As expected, for weak spin-orbit 
coupling (dashed line), the spin polarization grad- 
ually decreases into the interior of the stripe and 
vanishes completely at the center. A completely 
other behavior exhibits the solid line calculated for 
£1 = 1. A spin-coherent standing wave appears 
for strong spin-orbit coupling. The period of this 
pattern is comparable with the mean-free path / 




i = 2Ky 



FIG. 2: Spin polarization vf z as a function of (Q = Kl = 1). Such rapid oscillations are in prin- 
£ = 2Ky for = 0.1 (dashed line) and Q = 1 

(solid line). ciple beyond the accessibility of the drift-diffusion 

approach. However, as already mentioned above, 
similar oscillations appear also in the ballistic regime. Moreover, recently it has been stated 
that the macroscopic transport equations work well even when the spin-diffusion length is 
comparable to the mean-free path. [31 ] 

In the strong-coupling limit Q ^> 1, the solution in Eq. (|39|) simplifies to 

f z0 cos(4^/v / 6) 



f z (k,y,s) 



(41) 



s cos(4Q£o/ V6) ' 

The period of this oscillation [given by lyEir/ \A(Kl) 2 )] is much smaller than the mean-free 
path I. In addition, the spin polarization f z in Eq. (j4ip diverges, when the half width L of 
the stripe satisfies the resonance condition (2n + l)vr/2 = 4fi£ /v / 6 with n being any integer. 
The appearance of these divergencies indicates that elastic scattering alone cannot always 
balance the tendency of mutual spin alignment in the strong spin-orbit coupling regime. To 
avoid this defect of the approach, obviously both elastic and inelastic scattering have to be 
treated on a full microscopic level. 
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To give a rough idea of the temporal evolution of spin-coherent standing waves, we treat 

Eq. fl39|) in the quasi-stationary regime (sr <C 1, 
i.e., o"o = 1) and perform the inverse Laplace trans- 
formation by a numerical integration. Figure 3 
shows an example for f z clS du function of £ and 
t/r. Initially, at t = 0, there is only the provided 
spin polarization f zQ at the boundaries. The time 
evolution starts with an intricate dependence that 
is by no means correctly reproduced by the drift- 
diffusion approach because it appears at a time 
scale t < t. After this short period of rapid spatio- 
temporal changes, the periodic spin pattern grad- 
ually emerges. Whereas the gross features of this picture are certainly reliable, a detailed 
analysis of the time evolution requires an approach based on a microscopic theory. 




£ = 2Ky 



FIG. 3: Spin polarization f z as a function of 
t/r and £ = 2Ky for Q, = 2 and f z0 = 0.2. 



V. SUMMARY 



Starting from microscopic quantum-kinetic equations for the spin-density matrix, we de- 
rived drift-diffusion equations for a 2DEG with Rashba spin-orbit interaction. The approach 
is valid in the weak as well as strong coupling regime. In the absence of external electric 
and magnetic fields, the set of equations decouples. The first set of equations relates the 
charge density to the transverse in-plane spin component. Another set of equations links 
the out-of plane spin component with the longitudinal in-plane spin contribution. Whereas 
both sets of equations have an universal character in the weak scattering limit, they exhibit 
interesting peculiarities for strong spin-orbit coupling. In the latter case, the spin relaxation 
time t s becomes shorter than the elastic scattering time r. Consequently, spin effects have 
a principal non-Markovian character in this regime. Moreover, the diffusion coefficient D z 
of the f z -fd spin channel changes its sign with increasing spin-orbit coupling Q = Kl. This 
dependence gives rise to an instability in the spin subsystem (D z = 0) and to an undamped 
oscillatory spin pattern (D z < 0). The wavelength of the long-lived coherent spin rotation 
is smaller than the mean-free path. This fact conflicts with basic assumptions of the drift- 
diffusion approach, which is only applicable for times much longer than the elastic scattering 
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time and for diffusion lengths much larger than the mean-free path. The approach signals 

y studied transition from a uniform to 
321 ] The rapid spatial variations of the 



the existence of an instability similar to the recent 
a nonuniform ground state in the Rashba model, 
out-of plane spin polarization gives rise to microscopic circulating currents via the induced 
magnetic field. Within a self-consistent schema, one expects that the retroaction of these 
currents on spins may lead to a stabilization of the spin subsystem and to a finite damping 
of spin oscillations. 

In order to appreciate the unusual results obtained from the drift-diffusion equations in 
the strong-coupling limit, an exact treatment of the original quantum- kinetic equations is 
desirable. For strong spin-orbit coupling (fi > 1), the spin relaxation time is much smaller 
than the scattering time. Therefore, most interesting is the consideration of the ballistic 
spin regime. Exact analytical solutions of the kinetic equations confirm the existence of a 
long-lived oscillating spin pattern. This confirmation strongly suggests that spin-coherent 
standing waves in a semiconducting stripe are indeed serious solutions of the drift-diffusion 
equations. Moreover, a recent alternative approach based on Monte-carlo simulations also 
arrived at the conclusion that spin-coherent standing waves have an extremely long spin- 
relaxation time. {28] 

We conclude that the prediction of an oscillatory spin pattern in semiconductors with 
strong spin-orbit coupling of the Rashba type seems to be interesting both from a theoretical 
and experimental point of view. Further progress is expected from a due treatment of the 
quantum-kinetic equations for the density matrix beyond the pure ballistic regime (r — > 00). 
In addition, external electric and magnetic fields sensitively influence spin coherent waves by 
mixing all components of the density matrix. The experimental verification of the predicted 
long-lived spin-coherent waves could be useful for future spintronic applications. 
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